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ABSTRACT 

The fluxes of heavy nuclei with Z > 6 in the primary cosmic radiation were 
measured with a nuclear emulsion stack exposed in a balloon flight over Texas 
on March 26-27, 1962, and the fluxes of M(6 1 Z < 9) and H(Z > 10) nuclei extra- 
polated to the top of the atmosphere were found to be 6.20 ± 0.46 and 2.02 ± 0.25 
particles/m 2 ster sec, respectively. Results of this work and those of other in- 
vestigations made during the last solar cycle (1954-1964) are used to study the 
long term solar modulation of primary S-nuclei (Z > 6) of rigidity >4,8 BV. It is 
found that the lag between the primary intensity changes and the sunspot number 
is about one and a half years during the increasing phase of solar activity and is 
very small during the declining phase of solar activity. This is in agreement 
with the changes of the spectral shape during the solar cycle as suggested by 
Simpson. The implication of the time lag in terms of the scale size of the tran- 
sition region of the magnetic field responsible for the long term modulation is 
briefly discussed. 

*NASA— National Academy Postdoctoral Research Associate. Present address: Tata Institute of 
Fundamental Research, Bombay 5. 
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THE FLUX OF HEAVY NUCLEI IN THE 


PRIMARY COSMIC RAYS OVER TEXAS 
IN MARCH 1962, AND THE SOLAR CYCLE MODULATION 

INTRODUCTION 

The study of fluxes of heavy nuclei in the primary cosmic rays with energies 
£1.6 Bev/nucleon incident on Texas (geomagnetic latitude k % 41°) is of interest 
from two considerations. Firstly, from the comparison of the fluxes measured 
over Texas with those measured simultaneously at a low geomagnetic latitude, 
e.g. near the equator, one can obtain a two point integral energy spectrum of the 
different components and compare them in the energy range of about 1.6 to 7 
Bev/nucleon. Secondly, from the comparison of the fluxes of heavy nuclei meas- 
ured over Texas at different times over the solar cycle, the long term solar 
modulation of the primary heavy nuclei of rigidity £4.8 BV may be investigated. 
With these objectives in mind a nuclear emulsion stack was exposed over Texas 
on March 26-27, 1962. Simultaneously another emulsion stack was exposed over 
Hyderabad, India by the group at the Tata Institute of Fundamental Research, 
Bombay. The preliminary results from the Texas stack were reported previously 
(Neelakantan and Biswas, 1963). Using the initial results from the Texas stack 
and the Hyderabad stack, the integral energy spectrum of the heavy nuclei in the 
energy interval 1.6 to 7 Bev/nucleon was obtained (Badhwar et al. 1963). This 
aspect of the problem will not be discussed further in this paper. Here we pre- 
sent the final results of the analysis of the heavy nuclei in the Texas Stack, the 
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comparison of these flux values with those of other investigators obtained during 
the solar cycle (1954-64), and a discussion of the long term solar modulation of 
primary heavy nuclei of rigidity ^4.8 BV. 

EXPERIMENTAL PROCEDURE 

The experiment was conducted with an emulsion stack flown in a balloon 
flight over Texas on March 26-27, 1962. The balloon floated for 14.5 hours at 
5.55 g/cm 2 of residual atmosphere. The time-altitude curve and the trajectory 
of the balloon are shown in Figures la and lb, respectively. An arrangement to 
provide rotation of the stack by 90° upon reaching ceiling failed so that the flip- 
ping occured at the time of launch. Standard nuclear emulsion techniques were 
used in the analysis of the stack, the details of which are given in the appendix. 

Since the cutoff rigidity allows the tracks to be considered as relativistic, 
charge identification was made by at least two independent ionization measure- 
ments on each particle track. The ionization measurements were made by S-ray 
counting in the G5 emulsions and blob-gap counting in the G2 emulsions by three 
different observers. Using this technique to obtain an average Z- value, 185 M 
nuclei, and 66 H nuclei were detected. Correcting for scanning efficiency this 
gives 195.7 M nuclei and, 66 H nuclei at the scan line. 

RESULTS 

The observed fluxes of M and H nuclei were extrapolated to the top of the 
atmosphere by the following two methods. Firstly in the extrapolation procedure 
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using diffusion equations, we use the values of interaction mean free paths in 
air for M and H nuclei as 27.1 and 18.9 g/cm 2 respectively and the best available 
values of fragmentation parameters, P MM , P HM , and P HM obtained as 0.15, 0.32, 
and 0.23 respectively from the weighted average values measured in graphite, 
celluloid and teflon (Durgaprasad, 1963; Friedlander et al., 1963). These values 
are in general agreement with those used by other investigators (e.g., O'Dell 
et al., 1961) for an air like medium. The fluxes thus obtained after correcting 
for particles recorded during the ascent and extrapolated to the top of the at- 
mosphere are 

J M (0) - 6.20 ± 0.46 particles/m 2 ster sec 

J H (0) = 2.30 ± 0.28 particles/m 2 ster sec 

J s (0) = 8.50 ± 0.54 particles/m 2 ster sec 

It should be pointed out here that in case of H nuclei, the diffusion method 

predicts the absorption mean free path of H nuclei in air as A H = 24.5 g/cm 2 
which is smaller than that of M nuclei A M = 31.9 g/cm 2 , so that, the ratio of 
H/M in the atmosphere should increase with atmospheric depth. However, the 
observed H/M ratio as a function of the atmospheric depth is found to be slowing 
decreasing or almost constant as shown by Daniel and Durgaprasad (1962) and 
Webber (1965). As pointed out by these authors, the diffusion procedure probably 
breaks down in the case of H nuclei because of the fact that the group is so large. 
Therefore, the second procedure is used to obtain the extrapolated flux for H 
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nuclei. Thus the observed ratio of H/M at the flight altitude (after correcting 
for ascent particles) and the observed growth curve of H/M in air as given by 
Webber (1965) are used to obtain the (H/M) 0 ratio at the top of the atmosphere. 
Then the J H (0) value at the top of the atmosphere is determined by using the 
J M (0) as obtained previously. The fluxes thus obtained are 

J M (0) = 6.20 ± 0.46 particles/m 2 ster sec 

J H (0) = 2.02 ± 0.25 particles/m 2 ster sec 

J s (0) = 8.22 ± 0.52 particles/m 2 ster sec 

In any case, since the measurements were made at small atmospheric depth 
(5.55 g/cm 2 ), the above two different extrapolation procedures yield results 
which are similar. 

COMPARISONS WITH OTHER INVESTIGATIONS 

The available data on the flux of (M + H) nuclei measured over Texas at 
atmospheric depth of <12 g/cm 2 by different investigators are given in Table I. 
Before comparing these results, the following two points should be considered. 
Firstly, different investigators have used different methods and parameters for 
the extrapolation of the measured fluxes, particularly for the H-nuclei. The 
absorption mean free path of H-nuclei in air, A H , used by different authors vary 
widely from 24 to 54 g/cm 2 . For M-nuclei, however, most of the authors have 
used nearly the same value of A M as ~32 g/cm 2 . In order to compare different 
measurements correctly, it is important to normalize them using the same 
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extrapolation procedure. Therefore, we have obtained the J M ' (0) and J H ' (0) 
fluxes of different investigators using the same method of extrapolation, namely 
the second procedure given in the previous Section. The normalized fluxes 
J s ' (0) are shown in column 9 of Table I. Secondly, most of the investigators of 
the Texas experiments assume the geomagnetic cut-off energy as 1.5 Bev/nucleon 
or cut-off rigidity as 4.5 BV, irrespective of the trajectory of the balloon. We 
have calculated the effective cut-off rigidities for our balloon trajectory as well 
as those of other investigators who have given their balloon trajectories, taking 
into account the variation of the cut-off rigidity with geographic location during 
the flight as given by Quenby and Wenk (1962) and assuming an exponent 1.5 for 
the integral energy spectrum. The effective geomagnetic cut-off rigidities are 
also shown in Table I. Since the average cut-off rigidity for these flights is about 
4.8 BV, we have normalized J s ' (0) flux values to the same cut-off rigidity of 
4.8 BV. The final normalized fluxes, J s " (0), are shown in the last column of 
Table I and these are used for further analysis. 

SOLAR CYCLE VARIATION 

Next we consider the variation of the flux of primary heavy nuclei of Z > 6 
of rigidity >4.8 BV during the last solar cycle, 1954-1964. So far, a number of 
investigations have been made on the solar cycle variation of the flux of low 
energy primary protons and helium nuclei of rigidity 0.9-3 BV (See, for example. 
Summaries by Webber, 1962, 1965; Freier and Waddington, 1965) and relatively 
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less has been discussed on the long term modulation of heavy nuclei of rigidity 
>4.8 BV. In studying the long term solar modulation of the heavy nuclei we have 
taken the flux of M + H nuclei as a whole so that statistical errors are reduced 
and the uncertainties if any, in the separation of M and H nuclei in different 
experiments do not enter. 

The available data on the normalized flux of (M + H) nuclei of rigidity 
>4.8 BV, as given in Table I, have been plotted against Mt. Washington neutron 
monitor counting rate, in Figure 2. It is seen that most of the flux measurements 
were made during periods close to either solar maximum, (Mt. Washington rate 
about 2000) or to solar minimum (Mt. Washington rate about 2400) and very few 
data are available for intermediate periods. The data point from the present 
work is the only one available for the period of intermediate neutron monitor 
counting rate. The variation shown in Figure 2 is different from that shown by 
Webber (1965) because in the latter the preliminary value of the flux obtained in 
the present work was used. Also we have neglected the flux value measured by 
McDonald and Webber (1962) on 20th March, 1956, which was included in Webber's 
Figure 13 because this measurement was made during a period of large Forbush 
decrease (McDonald, 1957). Since in some Forbush decrease events the rigidity 
dependence is different from that due to the 11 year modulation (Webber, 1962; 
Kane, Wincker and Arnoldy, 1965) it is important that Forbush decrease events 
should not be included in the study of the long term variation. The regression 
curve given by Webber (1965) would predict the flux of (M + H) nuclei at the 
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time of this investigation as 7.0 particles/m 2 ster sec, which is more than two 
standard derivations lower than the measured value. Although a nearly linear 
regression curve cannot be ruled out completely on the basis of the single data 
point, the nonlinear regression curve shown in Figure 2 is indicative of the fea- 
ture that nuclei of rigidity >4.8 BV seem to recover faster than the Mt. Wash- 
ington neutron monitor which responds to all particles of rigidity greater than 
1 BV. This is consistent with the neutron monitor data available at present, as 
discussed below. 

In order to study in more details the changes of the heavy primary flux of 
rigidity >4.8 BV with the increasing and decreasing phases of the solar activity 
we have plotted in Figure 3, the fluxes of the S nuclei for the years 1954-1963, 
together with the solar activity in sunspot number in the inverse scale. As seen 
in the figure, during the increasing phase of the solar activity from 1954 to 
1957-58, a number of measurements of the S nuclei flux are available which show 
that between solar minimum in 1954 and solar maximum in 1957, the S nuclei 
flux decreased by as much as 45 percent below the 1954 level. It also shows 
that the decrease in the primary flux lagged behind the increase of solar activity 
(as given by the sunspot number), the amount of lag being a function of time. 
During the low and intermediate period of solar activity this lag was about one 
and one-half years and this lag seems to decrease as the maximum of the solar 
activity is reached. This behavior is of similar pattern to that shown earlier by 
Simpson (1963) and by others (Webber, 1965; Waddington, 1965) of the counting 
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rate of the Mt. Washington neutron monitor (responding to particles >1 BV). 

The lag of the latter with respect to sunspot number during the same period of 
the solar cycle was about one year which is nearly the same as that observed 
here for the S-nuclei. The minimum of the Mt. Washington neutron monitor and 
of the S nuclei flux both seem to occur very close to the maximum of solar 
activity. 

During the decreasing phase of the solar activity, there are unfortunately 
very few data on the flux of S nuclei of rigidity >4.8 BV, the present work being 
the only measurement available at present between the years 1958 to 1962. 

Hence, subject to this limitation, the inferences would be only suggestive and 
much in need of further verification. The fluxes measured in this experiment 
in 1962 and by Webber (1965) in 1963 are nearly the same within statistical errors 
indicating that by 1962-63 the flux of S nuclei recovered to a large extent and 
was only a few percent below the 1954 level. This suggests that during the de- 
creasing period of solar activity in 1958 to 1963, the flux of S nuclei seems to 
follow more closely the solar activity as given by the sunspot number, as opposed 
to the behavior of the high latitude neutron monitor which shows a large lag of 
about two years during 1960-1962. This is consistent with the general picture 
given by Simpson (1963) which indicates that during the recovery phase, neutron 
intensity changes produced by primary particles with rigidity greater than 2.8 
BV recovered faster than those produced by particles with rigidity greater than 
1.5 BV and the equatorial neutron counting rate produced by particles of rigidity 
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greater than 14 BV recovered almost fully to 1954 value during the end of 1961 
and early 1962. 

Another way of examining the lag is to plot the flux or its percentage de- 
crease below 1954 level against the mean sunspot number as shown in Figure 1. 
This shows the characteristic "hysteresis" effect (Simpson, 1963), and how the 
phase lag is changing over the solar cycle. Here again the lack of data points 
between the years 1958-1962 makes the recovery phase during the years 1958-62 
uncertain. The rigidity dependence of the modulation as discussed earlier would 
produce a hysteresis loop such that, the higher the rigidity, the smaller would be 
the area under the hysteresis loop. 

In the last sunspot cycle the onset of modulation for S nuclei of rigidity 
>4.8 BV appears to occur at a solar activity level given by a sunspot number of 
about 100, whereas the recovery for the same particles appears to have occurred 
by the time the sunspot number reached 50. Within the limits of the data these 
results agree with the estimates one would obtain from Simpson’s data (1963). 

Similar behavior can be seen in the a-particle fluxes measured over Texas 
during the solar cycle 1954-1963 as given by Freier and Waddington (1965). Here 
also the scatter of the points are large as statistical errors are considerable and 
there is a lack of data points between 1959 and 1963 so that the declining cycle is 
not well known. We wish to point out that from solar minimum (1954) to solar 
maximum (1957), the flux of a-particles of rigidity £4.8 BV decreased from 90 to 
65 particles/m 2 ster sec, i.e. a decrease of about 28% below 1954 level in contrast 
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to 45% for S nuclei. This difference is not understood at the present if it cannot 
be attributed to statistical fluctuations, as He and S nuclei have about the same 
mass to charge ratio and hence are expected to be modulated in a similar manner. 
Further studies are needed to clarify this effect. 

Thus we have the overall picture which seems to indicate that during the 
increasing phase of the solar activity, all primaries of rigidity >1 BV and >4.8 BV 
decrease in a somewhat similar manner with a lag of one to two years with sun- 
spot number, the lag being greater for particles of higher rigidity, while during 
the declining phase of solar activity, the corresponding lag decreases as the rigid- 
ity of the particle increases. This is consistent with the changes in the spectral 
shape of the primary radiation for the periods 1954-1958 and for 1958- 
1962 as given by Simpson (1963). The study of the changes of the spectral 
shape and the intensities over the entire solar cycle provides important 
clues regarding the origin and structure of the magnetic fields which lie 
beyond the orbit of the earth and one can estimate the scale size of the 
transition region from the characteristic build-up and decay times as dis- 
cussed by Simpson (1963). Unfortunately the data on S-Nuclei fluxes during 
the declining phase of solar activity are much too sparse to draw any con- 
clusions regarding the depth of the modulating region although the data are 
not inconsistent with the estimates made by Simpson based on neutron moni- 
tor studies over a range of rigidities. The estimates by Simpson and others by 
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Parker (1963) and Axford et al. (1963) based on termination of solar plasma 
by galactic magnetic field range from 10-90 A. U. 

Appendix 

The emulsion stack consisted of 69 Ilford G5 and 7 G2 emulsions of size 
20 x 15 x 0.06 cm. Each G5 emulsion was interleaved after every 10th G5 
emulsion in the stack. Each emulsion sheet was scanned on a line 5 cm away 
from the top edge of the emulsion for all tracks crossing the scan line and 
satisfying the following criteria: 

Ionization >16 times minimum 
Projected length in emulsion >1.9 mm 
Zenith angle 150°. 

Only tracks with projected length >2 mm and with zenith angle <45 degrees were 
finally selected for measurement. The primary tracks with Z > 5 were separated 
from slow particle tracks by following the tracks through several emulsions. 

Then the selected tracks were followed back to the top of the stack and those 
which came from interactions above the scan line were rejected. 

Since the geomagnetic cut-off energy over Texas is about 1.6 Bev/nucleon, 
all particles entering the stack can be considered relativistic. Hence their 
charge was determined by making two separate ionization measurements: (i) The 
long S-rays which projected at a distance >l/x from either side of the stack were 
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counted. About 300 S-rays were counted on each track to determine its charge. 
The S-ray counting criteria was kept constant by recounting on some standard 
tracks frequently, (ii) Blob-gap counting was made on all tracks which passed 
through a G2 emulsion and belonged to L(3 Z < 5) and M(6 £ Z < 9) groups of 
nuclei as determined from S-ray density measurements. About 300 blobs and 
300 gaps of length >1.5/x were counted on each track. Plate to plate normalization 
of the blob density in G2 plates were made in the usual manner. This method was 
adopted for an independent charge determination for these nuclei in order to 
obtain reliable separation of L and M groups of nuclei. The charge calibration 
was done in the standard manner by means of tracks showing charge indicating 
interactions. 

On all heavy nucleus tracks at least two independent charge determinations 
were made from measurements of at least two different observers. Thus out of 
286 tracks of heavy nuclei measured, four independent charge determinations 
were made on 29 tracks, three independent determinations on 135 tracks, and 
two independent determinations on 122 tracks. In general there was good agree- 
ment between these charge assignments, excepting a few cases. In these cases 
mean charge was assigned on the basis of the average of the independent de- 
terminations. The cross plot of the charges determined by S-ray density and 
blob gap method is shown in Figure 5. Thus using all the available measurements 
good charge resolution was obtained between L and M nuclei as well as between 
M and H nuclei (Z > 10). No attempts were made to determine individual charges 
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of nuclei belonging to the H group (Z 2 10). The charge distribution as obtained 
from all the measurements is shown in Figure 6. 

A total of 34 L nuclei, 185 M nuclei and 66 H nuclei were obtained in the 
scan length of 240.5 cm and in scanned area of 13.60 cm 2 . 

The scanning efficiencies of the five scanners who scanned the emulsions 
were determined by rescanning of the plates by a different scanner, and these 
were found to be 100%, 100%, 88%, 96% and 80% for M nuclei and 100% for H 
nuclei. No attempt was made to analyze the L nuclei and the scan criteria were 
not chosen to efficiently accept this group. The weighted average of the scanning 
efficiency for M nuclei was 94%. After correcting for scanning efficiency, the 
number of M nuclei was 195.7 and of H nuclei 66. 
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Instrumental cut-off 
Trajectory not given by authors 



Figure Captions 


Figure 1: (a) The time altitude curve of the balloon flight, (b) The balloon 
trajectory. 

Figure 2: A plot of the data in Table I giving the flux of S-Nuclei versus the 

Mt. Washington Neutron Monitor during the most recent solar cycle 


(1954-1964). 


Figure 3: 


Figure 4: 


Figure 5: 


Figure 6: 


The S-Nuclei flux (©) and the smoothed sunspot number (x) plotted as 
a function of time during the solar cycle. 

The flux of S-Nuclei plotted as a function of the smoothed sunspot 
number during the last solar cycle (1954-1964). 

A comparison of the charge assignments made on the basis of blob-gap 
counts versus those made from delta-ray measurements (N s ) 

The charge distribution obtained using the averaged z-measurements 
described in the text. No attempt was made to efficiently detect and 
analyze the charges less than 6. 
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